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^^ ■ Let p be a homogeneous polynomial of degree n in n variables, 

p{zi, . . . , Zn) = p{Z), Z G C". We call such a polynomial p H-Stable 
li p{zi, . . . , Zn) ^ O provided the real parts Re{zi) > 0,1 < i < n. 
t^ ■ This notion from Control Theory is closely related to the notion of 

Hyperbolicity used intensively in the PDE theory. 

The main theorem in this paper states that if p{xi, . . . ,Xn) is a 
^ ' homogeneous H-Stable polynomial of degree n with nonnegative co- 

J^ . efïicients; degp{i) is the maximum degree of the variable Xj, Cj = 

uiin(degp{i) , i) and 

Capip) = inf Pi^'^^-'-'^n) 

Xi>0,l<i<n Xl ■ ■ ■ Xn 

then the following inequality holds 

Qn /(J. — w Ci-l 

p{0,...,0)>Cap{p) n 



dxi . . . dxn ' ' ^zr.t- V C, 

^ ^'' 2<i<n ^ ' 



This inequality is a vast (and unifying) generalization of the Van 
der Waerden conjecture on the permanents of doubly stochastic ma- 
trices as well as the Schrijver- Vahant conjecture on the number of 
perfect matchings in /c-regular bipartite graphs. These two famous 
results correspond to the H-Stable polynomials which are products 
of hnear forms. 

Our proof is relatively simple and "noncomputational" ; it uses just 
very basic properties of complex numbers and the AM/GM inequahty. 

1 The permanent, the mixed discriminant, 
the Van Der Waerden conjecture (s) and 
homogeneous polynomials 

Recall that an n x n matrix A is called doubly stochastic if it is nonnegative 
entry-wise and its every column and row sum to one. The set oi nxn doubly 
stochastic matrices is denoted by fi„. Let A{k,n) denote the set oi n x n 
matrices with nonnegative integer entries and row and column sums all equal 
to k. We define the foUowing subset of rational doubly stochastic matrices: 
ilk^n = {k~'^A : A e A{k, n)}. In a 1989 paper [2J R.B. Bapat defined the set 
Dn of doubly stochastic n-tuples of n x n matrices. 

An n-tuple A = {Ai, . . . , A„) belongs to Dn iff A >z O, i.e. Ai is a positive 
semi-definite matrix, 1 < i < n; trAi = 1 for 1 < i < n; ^27=1 ^i = -^? where 
/, as usual, stands for the identity matrix. Recall that the permanent of a 
square matrix A is defined by 

n 

per{A)= J2 UMh^ii))- 

Let us consider an n-tuple A = (Ai, j42, . . . A„), where Ai = {Ai{k,l) : 1 < 
A;, / < n) is a complex nxn matrix {l < i < n). Then 

DetA(ti,...,t„) = det( Y. ti^i) 

l<i<n 

is a homogeneous polynomial of degree n in ti,t2, ■■■ ,tn. The number 

D{A) := D{Au A2,...,An) = j-^^Detj,{0, . . . , 0) (1) 

Oti ■ ■ ■ Oln 



is called the mixed discriminant oiAi,A2,..., A^. 

The mixed discriminant is just another name, introduced by A.D. Alexan- 
drov, for 3-dimensional Pascal's hyperdeterminant. The permanent is a par- 
ticular (diagonal) case of the mixed discriminant. I.e. define the foUowing 
homogeneous polynomial 

Froö?A(ti,...,t„)= n E M^^J)h- (2) 

l<i<n l<j<n 

Then the next identity holds: 

gn 

Per{A) = ProdAiO, . . . , 0). (3) 

Oti , . . . , Otn 

Let us recall two famous results and one recent result by the author. 



1. Van der Waerden Conjecture 

The famous Van der Waerden Conjecture [23] states that 

min per (A) = — =: vdw{n) (VDW-bound) 

AGiin Tl 

and the minimum is attained uniquely at the matrix J„ in which every 
entry equals K The Van der Waerden Conjecture was posed in 1926 
and proved in 1981: D.I. Fahkman proved in [5j the lower bound ^; the 
fuU conjecture, i.e. the uniqueness part, was proved by G.P. Egorychev 
in H]. 

2. Schrijver- Valiant Conjecture 

Define 

X{k,n) = mm{per{A) : A G flk,n} = fc~" min{per(A) : A E A{k,n)}; 

eik) = hm(A(A;,n))". 

n^oo 

It was proved in [26] that, using our notations, 9{k) < G{k) =: (^^)^^^ 
and conjectured that 9{k) = G{k). Though the case of A; = 3 was 
proved by M. Voorhoeve in 1979 [28], this conjecture was settled only 
in 1998 [27] (17 years after the pubhshed proof of the Van der Waerden 



Conjecture). The main result of [27| is the remarkable (Schrijver- 
bound): 

(k- ix^''-^)" 

mm{per{A) : A G f]fc,„} > ( —^ J (4) 

The proof of (Schrijver-bound) in [27j is, in the words of its author, 
"highly comphcated" . 

Remark 1.1: The dynamics of research which led to (Schrijver- 
bound) is quite fascinating. \ï k = 2 then mm^gA(2,n)Per(74) = 2. 
Erdos and Renyi conjectured in 1968 paper that 3-regular case already 
has exponential growth: 

min per(A) > a^,a > 1. 

AGA(3,n) 

This conjecture is imphed by (VDW-bound), this connection was 
another important motivation for the Van der Waerden Conjecture. 
The Erdos-Renyi conjecture was answered by M. Voorhoeve in 1979 



min per(A) > 6 - . (5) 

Amazingly, the Voorhoeve's bound ([5]) is asymptotically sharp and the 
proof of this fact is probabihstic. In 1981 paper [2ü], A. Schrijver and 
W.G.Vahant found a sequence ^k,n of probabihstic distributions on 
A{k,n) such that 

il^ (Ae^i.^n)^"'"^^)) " - ii^ {E,,,„per{A)) ~ = k l^-^-j (6) 

{LM. Wanless recently extended in ÏSÜf the upper bound ^ to the 
boolean matrices in A{k,n).) 



It follows from the Voorhoeve's bound ([5]) that 

hm [E,., per(A))" = hm ( min per(A)] for A; = 2,3. 

This was the rather bald intuition that gave rise to the Schrijver- Valiant 
1981 conjecture. 



/ r._, \ fc— 1 

The number ki—^j in Schrijver- Valiant conjecture came up via 
combinatorics foliowed by the Standard Stirling's formula manipula- 
tions. On the other hand G{k) = (^)'=-^ = ^j^^- 
I 

3. Bapat's Conjecture (Van der Waerden Conjecture for mixed 
discriminants) 

One of the problems posed in [2] is to determine the minimum of mixed 

discriminants of doubly stochastic tuples: minAeD„D{A) =? 

Quite naturally, R.V.Bapat conjectured that min^eD^ D{A) = ^ (Bapat- 

bound) and that it is attained uniquely at J„ =: (-/, . . . , -/). 

In [2] this conjecture was formulated for real matrices. The author 

proved it [13] even for the complex case, i.e. when matrices Ai above 

are complex positive semidefinite and, thus, hermitian. 

1.1 The Ultimate Unification (and Simplification) 

Falikman/Egorychev proofs of the Van Der Waerden conjecture as well our 
proof of Bapat's conjecture are based on the Alexandrov inequalities for 
mixed discriminants [T] and some optimization theory, which is rather ad- 
vanced in the case of the Bapat's conjecture. They all rely heavily on the 
matrix structure and essentially of non-inductive nature. 
(D. I. Falikman independently rediscovered in f^ the diagonal case of the 
Alexandrov inequalities and used a clever penalty functional. The very short 
paper 15^ is supremely original, it cites only three references and uses none 
of them.) 

The Schrijver's proof has nothing in common with these analytic proofs; it 
is based on the finely tuned combinatorial arguments and multi-level induc- 
tion. It heavily relies on the fact that the entries of matrices A G A(/c, n) are 
integers. 

The main result of this paper is one, easily stated and proved by easy 
induction, theorem which unifies, generalizes and, in the case of (Schrijver- 
bound), improves the results described above. This theorem is formulated 
in terms of the mixed derivative „ , ^" p(0, . . . , 0) (rewind to the formula 
Q) of H-Stable (or positive hyperbolic) homogeneous polynomials p. 
The next two completely self-contained sections introducé the basics of sta- 
ble homogeneous polynomials and proofs of the theorem and its coroUaries. 



We have tried to simplify everything to the undergraduate level, inaking the 
paper longer than a dry technical note of 4-5 pages. Our proof of the unique- 
ness in the generahzed Van der Waerden Conjecture is a bit more involved, 
as it uses Garding's resuh on the convexity of the hyperbohc cone. 

2 Homogeneous Polynomials 

The next definition introduces key notations and notions. 
Definition 2.1: 

1. The hnear space of homogeneous polynomials with real (complex) coef- 
ficients of degree n and in m variables is denoted HomR{m, n) {Homc{m, n)). 
We denote as Hom^{m,n) {Hom^j^{n,m)) the closed convex cone of 
polynomials p G HomR{m, n) with nonnegative (positive) coefïicients. 

2. For a polynomial p G Hom^{n,n) we define its Capacity as 

Cap{p) = inf p{xi, . . . , x„) = inf J'' " ' ' " ■ (7) 

3. Consider a polynomial p G Homc{m, n), 

P[Xi,...,Xm) = / , öri,...,rm _[J_ ^Jj . 
{ri,...,rm) l<i<m 

We define Rankp{S) as the maximal joint degree attained on the subset 
Sc{l,...,m}: 

Rankp{S) = max X!^i- (^) 

If S* = {i} is a singleton, we define degpii) = Rankp{S). 

4. Let p E Hom^{n,n), 

P[Xi, . . . , Xn) = / ^ CLri,...,r„ J_J_ ^i ■ 

riH \-r„=n l<i<n 

Such a homogeneous polynomial p with nonnegative coefficients is called 
doubly-stochastic if 

— p(l,l,...,l) = l:l<ï<n. 



In other words, p G Hom^{n, n) is doubly-stochastic if 

Y. «r-i,...,r„^j = 1 : 1 < i < n. (9) 



riH hr„=n 



It foUows from the Euler's identity that p(l,l,...,l) = 1: 

E «n,...,r„ = 1 (10) 

r-iH hr„=n 

Using the concavity of the logarithm on R++ we get that 
log(p(xi,...,x„)) > ^ a^i,...,r„ E ri\og{xi)=\og{xi---Xn). 

Therefore 

Fact 2.2: If p G Hom+{n,n) is doubly-stochastic then Cap{p) = 1. 

5. A polynomial jo G Homc{m, n) is called H-Stable iip{Z) ^ O provided 
i?e(Z) > 0; is called H-SStable if p{Z) ^ O provided Re{Z) > O and 
J2i<i<m Reizi) > 0. 

W^e coined the term 'H-Stable" to stress two things: Homogeniety and 
Hurwitz' stability. Other terms are used in the same context: Wide 
Sense Stable in 115], Half-PIane Property in |^. 

6. We define 

vdw{i) = %G(^ = -^^p^, = C-^T\^ > i;Gii) = 1. (11) 

Notice that vdw{i) as well as G{i) are strictly decreasing sequences. 
I 
Example 2.3: 

1. Letpe Hom+{2,2),p{xi,X2) = jxl + CxiX2 + ^xl]A,B,C>0. Then 
Capip) = C + y/AB and the polynomial p is H-Stable iff C > VAB. 

2. Let v4 G fin be a doubly stochastic matrix. Then the polynomial ProdA 
is doubly-stochastic. Therefore Cap{ProdA) = 1. In the same way, if 
A G Dn is a doubly stochastic ri-tuple then the polynomial DetA is 
doubly-stochastic and Cap{Detj^) = 1. 



3. Let A = {Al, A2, . . . Am) be an m-tuple of PSD hermitian n x n ma- 
trices, and J2i<i<m^i ^ O (the sum is positive-definite). Then the 
determinantal polynomial DetAih, . . . ,tm) = det(X]i<j<m^iA) is H- 
Stable and 

Rankoet^ {S) = Rank{J2 A) • (12) 

I 

The main result in this paper is the following Theorem. 

Theorem 2.4: Let p G Homj^{n,n) be H-Stable polynomial. Then the 
following inequality holds 

r. ^" p(0,...,0)> n G{min{i,degp{i)))Cap{p). (13) 

OXi . . . OXn 2<i<n 

Note that 

Y[ G{m.m{i,degp{i))) > Y\_ G{i) = vdw{n), 

2<i<n 2<i<n 

which gives the next generahzed Van Der Waerden Inequality: 

Corollary 2.5: Let p E Homj^{n,n) be HStahle polynomial. Then 

^^—^j(,....,,)>-Ca,(p). (14) 

Corollary ^2. 51) was conjectured by the author in fTD^ . where it was proved 

that 

dx ^"dx P^^' • • • 5 0) — C{n)Cap{p) for some constant C{n). 



2.1 Three Conjectures/Inequalities 

The fundament al nature of Theorem (12.41) is illustrated in the following Ex- 
ample. 

Example 2.6: 



1. Let j4 G r2„ be n X n doubly stochastic matrix. It is easy to show that 
the polynomial PtocLa is H-Stable and doubly-stochastic. Therefore 
Cap{ProdA) = 1- Applying Corollary (12.51) we get the celebrated Fa- 
hkman's result |5]: 



min per(A) = — . 



( The complementary uniqueness statement for Corollary 1^2. 5\} will be 
considered in Section^.) 



2. Let [Al, . . . , An) = A G Dn be a doubly stochastic ra-tuple. Then the 
determinantal polynomial Det^ is H-Stable and doubly-stochastic. 
Thus Cap{Detfi^ = 1 and we get the (Bapat-bound), proved by the 

author: 

n\ 
min -D(A) = — . 

AeDn ^ ' n" 

3. Important for what foUows is the next observation, which is a diagonal 
case of (fT2|): 

degprod^ij) is equal to the number of nonzero entries in the jth 
column of the matrix A. 

The next Corrolary combines this observation with Theorem (l2.4l) . 

Corollary 2.7: 

(a) Let Cj be the number of nonzero entries in the jth column of A, 
where A is an n x n matrix with non-negative real entries. Then 

per{A)> n G{mm{j,Cj))Cap{ProdA). (15) 

2<i<n 

(b) Suppose that Cj < k : k + 1 < j < n. Then 

per{A) > ( (-^) j -^CapiProdA). (16) 

Let A(/c, n) denote the set of n x n matrices with nonnegative integer 
entries and row and column sums all equal to k. The matrices in A{k, n) 



correspond to the A;-regular bipartite graphs with multiple edges. 
Recall the (Schrijver-bound): 



min per{A) > k'^G^kY 

AeA{k,n) 



{k-iy 



k-l\ " 



k 



fc-2 



(17) 



The Falikman's inequality gives that 

min per{A) > k'^vdw{n) > k'^G{kY if k>n. 

AeA{k,n) 



Therefore the inequality fITTj) is interesting only ii k < n. 

Note that if A G A(A;, n),k < n then all columns of A have at most k 

nonzero entries. 

If ^ G A(A;, n) then the matrix ^A G fin, thus Cap(ProdA) = A;". As we 

observed above, degprodAÜ) — ^- Applying the inequality (TTüj) to the 

polynomial ProdA we get for A; < n an improved (Schrijver-bound): 



min per{A) > A;" 

A£A{k,n) 



k- l\ k-i\ " k\ 



k^ 



k 



> 



\k-l) 



fc-i- 



k 



fc-2 



(18) 



Interestingly, the inequality (TTSl) recovers for A; = 3 the Voorhoeve's 
inequality ([5]). 

4. The inequality (fT5|) is sharp if Cj = . . . = C„_i = n; Cn = k : 1 < k < 
n — 1. To see this, consider the doubly stochastic matrix 



D 



( a 

a 
c 

V c 



a b \ 

a b 

c O 

c O J 



l-b _ k-1 1 _ 1 

ra — 1 A;(n — 1)' A;' n — l' 



(19) 



and the associated polynomial 



Produi^xi, . . . ,Xr. 



Y^ axi) + bxn ( X! ^^i)' 



n—k 



l<i<n-l 



l<i<n-l 



10 



Since the matrix D is doubly stochastic, Cap{Produ) = 1- Direct 
inspection shows that 

fr) — Tl' 
per{D) = {n- l)\{kh)a^'\''-^ = G{k)- 



(n- l)"-i 
Which gives the equahty 

per{D) = Cap{Prodi)) Yl ^ i^^Ui ^j)) ■ 

2<j<n 

It follows that min {per (A) : A e ü^^^ = (£ïj^ (^)" ^ ^here Vt^^^ 
is the set of n x n doubly stochastic matrices with at least one zero 
entry. 



I 



2.2 The Main Idea 

Let p G Hom+{n, n). Define the following polynomials qi G Hom+{i, ï): 

Qin—i 

Qn =p,qi{xi,...,Xi) = — p{xi,...,Xi,0,...,0);l <i<n-l. 

OXi^l . . . OXn 

Notice that qi{xi) = q^^^-q^p{^)^i aiid 



Qn 1 ( d"^ 9" 

dxi . . . dxn 2 \ (9a;i(9a;i . . . dxn ^ 8x28x2 ■ ■ ■ dxr, 

(20) 

^' axi...axn' 



2 



^'^"^''"^'^ " d^i — a^^^°^'^''^'^2 [f^^^ ^n-Pi^)^i + inniz KrPi^)^2 

Therefore, Cap{qi) = jf-rifTPiO) and 



dXi . . . dXn V ÖXiÖXi . . . dXn 8X28X2 ■ . . dXn 

(21) 
Define the univariate polynomial R{t) = p{xi, . . . ,Xn-i,t)- Then its 
derivative at zero is 

R'iO)=qn^iix,,...,Xn^i). (22) 

Another simple but important observation is the next inequality: 

degqii) < min {i,degp{i)) <^=^ G {degg.{i)) > G {m.m{i, degp{i))) : 1 < i < n. 

(23) 

11 



Recall that vdw{i) = %. Suppose that the next inequalities hold 

Cap{q,^,) > Cap{qi) f'l'.^'\, = Cap{q;)G{{) : 2 < i < n. (24) 
vaw[z — 1) 

Or better, the next stronger ones hold 

Capiq^.i) > Cap{qi)G {deg^M : 2 < ï < n, (25) 

where 

G{m)= ^ ^ = . 26 

vaw\m — 1) \ m / 

The next result, proved by the straigthforward induction, summarizes the 
main idea of our approach. 

Theorem 2.8: 

1. If the inequalities [24\ ) hold then the next generalized Van Der Waerden 
inequality holds: 

p{0, . . . , 0) = Cap{qi) > vdw{n)Cap{p). (27) 



dxi . . . dxn 
In the same way, the next inequality holds for Cap{q2): 



gn / Qn Qn 

Pi^)+\h^^ ^^P(^)^^^ ^P(O) > 2vdw{n)Cap{p). 



dXi . . . dXn V ÖXiÖXi . . . dXn 9X20X2 . . . ÖX. 

(28) 

2. Ifthe inequalities [2B^] hold then the next generalized (Schrijver-bound) 
holds: 

Qn 

^— P(0' • • • ' 0) = Cap{qi) > Cap{p) Y[ G{ min(i, degp{i))). 

OXi . . . OXn 2<i<n 

(29) 

What is left is to prove that the inequahties (l25l) hold for H-Stable 
polynomials. 
We break the proof of this statement in two steps. 



12 



1. Prove that if p G Hom^{n,n) is H-Stable then g„_i is either zero or 
H-Stable. Using equation (1221) . this implication follows from Gauss- 
Lukas Theorem. Gauss-Lukas Theorem states that ii zi, . . . ,Zn & C are 
the roots of an univariate polynomial Q then the roots of its derivative 
Q' belong to the convex huil CO{{zi, . . . , Zn}). 

This step is, up to minor perturbation arguments, known. See, for 
instance, JT^ . The result in JT^ is stated in terras of hyperbolic poly- 
nomials, see Remark li5.S\) for the connection between H-Stable and 
hyperbolic polynomials. Our treatment, described in Section^, is self- 
contained, short and elementary. 

2. Prove that Cap{qn-i) > G{degp{n))C ap{p) . This inequality boils down 
to the next inequality for the univariate polynomial R from (!22|) : 

i?'(0)>G(de^(i?))(^inf^^. 

We prove it using AM/GM inequality and the fact that the roots of 
the polynomial R are real. 

R is instructive to see what is going on in the "permanental case": we 
start with the polynomial Prod^ which is a product of nonnegative linear 
forms. The very first polynomial in the induction, g„_i, is not of this type 
in the generic case. I.e. there is no one matrix/graph associated with qn-i- 
We gave up the matrix structure but had won the game. 

In the rest of the paper Facts are statements which are quite simple and 
(most likely) known. We included them having in mind the undergraduate 
student reader. 



3 Univariate Polynomials 

Proposition 3.1: 

1. (Gauss-Lukas Theorem) 

Let R{z) = J2o<i<n ö^j-2* be a Hurwitz polynomial with complex coeffi- 
cients, i.e. all the roots of R have negative real parts. 
Then its derivative R' is Hurwitz. 

2. Let R{z) = J2o<i<n öj-2* be a Hurwitz polynomial with real coefficients 
and a„ > 0. Then all the coefficients are positive real numbers. 

13 



Proof: 

1. Recall that 



R'(Z) y 



Let /i be a root of K . Consider two cases. First: yU is a root of R. Then 
clearly Re{ii) < 0. Second: /i is not a root of R. Then 

Suppose that Re{fi) > 0. As {a + ib)~^ = -^-^ we get that 

\fJ. - zj {Re{fi) - Re{zj)f + {Im{^i) - Im{z^)f 

Therefore Re{L) > O which leads to a contradiction. Thus Re{fi) < O 
and the derivative R' is Hurwitz. 

2. This part is easy and well known. 



The next simple result binds together all the small pieces of our approach. 

Lemma 3.2: Let Q{t) = Y.o<i<k c^j^*! Uk > 0,k > 2 be a polynomial with non- 
negative coefficients and real (non-positive) roots. Define C = infj^o^^- 
Then the next inequlity holds: 

ai = Q'(0) > ( -^ j C. (30) 

The equality holds if and only if all the roots of Q are equal negative numhers, 
i. e. Q{t) = b(t + a)'^ for some a,b > 0. 

Proof: If g(0) = O then Q'(0) >C> {^^f'^ C. 

Let Q{0) > 0. We then can assume WLOG that Q{0) = 1. In this case all 

the roots of Q are negative real numbers. Thus 



Q{t) := Yliait + l):ai>0,l<i<k, 



i=l 

14 



and Q'(0) = ai + . . . + ttfc. 

Using the AM/GM inequality we get that 

Ct < Qit) < P{t) =: M + ^y^t ] , t > 0. (31) 

It follows from basic calculus that 

inf ^ = P{s) = g'(0) {j^\ , where s - ^ 



t>o t ' ' ^ ' ' V^-i; ' Q'{Q){k-l) 

Therefore 



C<inf:^ = g'(0) 



k-l 



t>o t ^ ' \k-lj 

which finally yields the desired inequahty 

Q'{0)>('^] C,k>2. 

It follows from the uniqueness condition in the AM/GM inequality that 
the equality in (J3Ü1) holds if and only if O < ai = . . . = Ofc. | 

Remark 3.3: The condition that the roots of Q are real can be relaxed in 

several ways. For instance the statement of Lemma (13.21) holds for any map 

ƒ : R+ -^ R+ such that the derivative /'(O) exists and ƒ ^ is concave. 

If such map is log-concave, i.e log(/) is concave, then /'(O) > ^ infi>o ^■ 

Notice that the right inequality in flïïTl) is essentially equivalent to the con- 

cavity 

of the function ((5(t))'= on i?+. 

It was shown in [TT] that the inequality fl5Dl) is equivalent to the (VDW- 

bound) for doubly stochastic matrices A E fin '■ ^ = [ö^I^I • • • I^] with two 

distinct columns. | 

4 Stable homogeneous polynomials 

4.1 Basics 

Definition 4.1: Apolynomialp G Horricim,!!,) is called H-Stable if p(Z) 7^ 
O provided Re{Z) > 0; is called H-SStable if p{Z) ^ O provided Re{Z) > O 
and Ei<i<m Reizi) > 0. | 

15 



Fact 4.2: Let p G Homdfn, n) be H-Stable and A is m x m matrix with 
nonnegative real entries without zero rows. Then the polynomial pa, defined 
as Pa{Z) = p{AZ) is also H-Stable. If all entries of A are positive then pA 
is H-SStable. 

Fact 4.3: Let p E Homc{m,n), Y G C"^,p{Y) ^ 0. Define the following 
univariate polynomial of degree n: 

Lx,Y{t) = p{tY -X)= piY) \{{t- K;y{X)) -.XeC^. 

l<i<n 

Then 

Xi.YibX + aY) = bXi.YiX) + a;piX)=p{Y) H Xi-riX). (32) 

l<j<n 

The following simple result substantially simplifies the proofs below. Propo- 
sition (14.41) connects the notion of H-Stability with the notion of Hyper- 
bolicity, see more on this connection in Subsection(5.1). 

Proposition 4.4: A polynomial p G Homdm, n) is H-Stable if and only 
if p{X) 7^ O : X G K^^ and the roots of univariate polynomials P{tX — Y) : 
X, y G i?™^ are real positive numhers. 

Proof: 

1. Suppose that p{X) 7^ O : X G R^_^ and the roots of univariate poly- 
nomials p{tX — y) : X, y G -R++ are real positive numbers. It fol- 
lows from identities fl32|) (shift L ^ L + aX > 0) that the roots of 
P{tX — L) : X G R'^_^,L G R"^ are real numbers. We want to prove 
that this property implies that p G Homc{m,n) is H-Stable. Let 
Z = Re{Z) + iIm{Z) G C^ : Im{Z) G /2™,0 < Re{Z) G /2™+. If 
p{Z) = O then also p{—iRe{Z) + Im{Z)) = O, which contradicts the 
real rootedness of p{tX - F) : X > O, F G i?" 



jm 



2. Suppose that p G Homc{rn,n) is H-Stable. Let X,Y E R^_^ and 

p{zX — Y) = O, 2; = a + &z. We need to prove that 6 = and a > 0. 
lih ^ Q then p{aX -Y + biX) = {hiYp{X - h-^i{aX - y)) ^ O as 
the real part Re{X — h^^i{aX — y)) = X > 0. Therefore b = 0. If 
a < O then -(aX - Y) E R'^^. Which implies that p{aX - Y) = 
(-l)"p(-(aX - Y)) ^ 0. Thus a > 0. 
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I 

We will use the following corollaries: 

Corollary 4.5: If Re{Z) G R^ and a polynomial p is H-Stable then 

\piZ)\>\piReiZ))\. (33) 

Proof: Since p is continuous on C™ hence it is sufficiënt to assume that 

Re{Z) e R"^^. 

It follows from identities (1321) that 



piZ) = p{ReiZ) + ilmiZ)) = p{Re{Z)) JJ (1 + ^A,), (34) 

i<j<'n. 

where (Ai,...,A„) are the roots of the univariate polynomial p{tRe{Z) — 
Im{Z)). Because Re{Z) G -R++, all these roots are real numbers. 

Therefore \p{Z)\ = |p(i?e(Z))| ni<,<„ |1 + ïA,| > \p{Re{Z))\. | 

Corollary 4.6: Let p G Homc{m,n) be H-Stable; X,Y e R"^ and O < 

X + y G -R++- Then all the roots of the univariate polynomial equat 
p{tX + Y) = are real numbers. 



'wn 



Proof: Let p{tX + F) = O, then also p{{t - l)X + {X + Y)) = 0. Since 
X + Y > O hence t — 1 ^ 0. As the polynomial p is homogeneous therefore 
p{X + (1 - ty^{X + Y)) = 0. It follows that (1 - t)~^ is real, thus t is also 
a real number. | 

Fact 4.7: Let p G Homc{m,n) be H-SStable (H-Stable). Then for all 
X G R^j^ the coefficients of the polynomial q = -^ are positive (nonnega- 
tive) real numbers. 

Proof: We prove ffist the case of H-SStable polynomials. 
Since q{X) = 1 we get from (132!) that q{Y) is a positive real number for 
all vectors Y G R^^. Therefore, by a Standard interpolation argument, the 
coefficients of q are real. We will prove by induction the following equivalent 
statement: if q E HorriR^m, n) is H-SStable and q{Y) > O for all Y G R^^^ 
then the coefficients of q are all positive. Write q{t; Z) = J2o<i<n't^Qi{Z)j 
where Z G C"""^, the polynomials qi G Hompi{m — l,n — 'i),0<ï<n — 1 
and qn{Z) is a real number. Let us fix the complex vector Z such that 
Re{Z) G R"^-^ and Re{Z) ^ 0. Since q is H-SStable hence all roots of 
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the univariate polynomial q{t; Z) have negative real parts. Therefore, using 
the first part of Proposition fl3.ip . we get that polynomials qi '■ O < i < n 
are all H-SStable. Since the degree of g is ??, hence qn{Z) is a constant, 
qn{Z) = g(l; 0) > 0. Using now the second part of Proposition fl3.ip . we see 
that qi{Y) > O for all Y G -R++ and O < i < n. Continuing this process 
we will end up with either m = 1 ot n = 1. Both those cases have positive 
coefïicients. 

Let p 6 Homc{m, n) be H-Stable and A> Q is m^m matrix with pos- 
itive entries such that AX = X. Then for all e > O the polynomials qi+eA G 
Homfi{'m,n), defined as in Fact fl4.2p . are H-SStable and limg^o^z+eyi = <?• 
Therefore the coefïicients of q are nonnegative real numbers. | 

From now on we will deal only with the polynomials with nonnegative 
coefïicients. 

Corollary 4.8: Let pi G Hom^{m,n) be a sequence o/H-Stable polynomi- 
als and p = lirai^oo Pi- Then p is either zero or H-Stable. 

Some readers might recognize Corollary Iji4-^ '^s a particular case of A. Hur- 
witz's theorem on limits of sequences of nowhere zero analytical functions. 
Our proof below is elementary. 

Proof: Suppose that p is not zero. Since pG ifom+(m, n) hence p(a;i, ... ,0;^) > 
O if Xj > O : 1 < j < m. As the polynomials pi are H-Stable there- 
fore \pi{Z)\ > \pi{Re{Z)) I : Re{Z) G R^^. Taking the limits we get that 
\p{Z)\ > \p{Re{Z)) I > O : Re{Z) G /?™+, which means that p is H-Stable. 
I 

Fact 4.9: For a polynomial p G Homc{m, n) we define a polynomial 
q G Horricirn — 1, n — 1) as 

d 

q{xi, . . . , x„_i) = - — p{xi, ..., Xm_i, 0). 

OXm 

Then the next two statements hold: 

1. Let p G Hom^{m,n) be H-SStable. Then the polynomial q is also 
H-SStable. 

2. Let p G Hom+{m,n) be H-Stable. Then the polynomial q is either 
zero or H-Stable. 



Proof: 

1. Let p G Hom+{m^n) be H-SStable and consider an univariate poly- 
nomial 

R[z) = piY- z):zeC,Y e C"'-\ 

Suppose that O 7^ Re{Y) > 0. It follows from the definition of H- 
SStability that R{z) 7^ O if Re{z) > 0. In other words, the univariate 
polynomial R is Hurwitz. It follows from Gauss-Lukas Theorem that 

qiY) = R'iO) ^ O, 

which means that q is H-SStable. 

2. Let p G Hom^{m,n) be H-Stable and q y^ 0. Take a.n m x m ma- 
trix A > 0. Then the polynomial pi+eA, Pi+tA{Z) = p{{I + eA)Z) 
is H-SStable for all e > 0. Therefore, using the first part, qi+eA is 
H-SStable. Clearly lim^^o Qi+eA = Q- Since g 7^ O, it follows from 
Corollary (gS]) that q is H-Stable. 

I 

Theorem 4.10: Let p E Hom+{n,n) be H-Stable, and 

d 
qn-i{xi, . . . , x„_i) = - — p{xi, ..., a;„_i, 0). 

ÖXn 

Then 

Cap{qn-^i) > Cap{p)G {degp{n)) . (35) 

Proof: We need to prove that 
d 



dXn 



p{xi,...,Xn-i,0)>Cap{p)G{degp{n)),Xi,...,Xn^i>0, Yl 2:^ = 1. 



l<i<n-l 



Fix a positive vector (xi, . . . , Xn-i), Y[i<i<n-i ^i = ^ and define, as in proof 

ofFactdMD, 

the polynomial R{t) = p{xi, . . . , a;„_i, t). It follows from Corollary fl4.6l) that 

all the roots of R are real. Since the coefïicients of the polynomial R are 

non-negative hence its roots are non-positive real numbers. It follows from a 

definition of Cap{p) that R{t) > Cap{p)t, therefore 

M^>Cap(p). 
t>o t ~ 
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The degree of the polynoinial R is equal to degp{n). It finally follows from 
Lemma fl3.2l) that 

Qn-iixi, ..., x„_i, 0) = R'{0) > Cap{p)G {degp{n)) . 



5 Uniqueness in Generalized Van Der Waer- 
den Inequality 

5.1 Hyperbolic Polynomials 

The following concept of hyperbolic polynomials arose from the theory of par- 
tial differential equations |^, JT^y . A recent paper ^241 gives nice and concise 
introduction to the area (with simplified proofs of the key theorems) and de- 
scribes connections to convex optimization. 

Definition 5.1: 

1. A homogeneous polynomial p : C™ —>■ C oi degree n [p & Homc{m., n)) 
is called hyperbohc in the direction e G R"^ (or e-hyperbohc) if p(e) ^ O 
and for each vector X G R"^ the univariate (in A) polynomial p{X — Xe) 
has exactly n real roots counting their multiplicities. 

2. Denote an ordered vector of roots of p{x — Ae) as 

Ae(X) = (A„(X) > A„_i(X) > . . . Ai(X)). 

Call X eR"" e-positive (e-nonnegative) if Ai(X) > O (A„(X) > 0). We 
denote the closed set of e-nonnegative vectors as Ne{p), and the open 
set of e-positive vectors as Ce{p). 



Remark 5.2: Proposition f l4.4p essentially says that a polynomial p G 
Homc{m, n) is H-Stable iff p is hyperbolic in some direction e G -R++ and 
the inclusion R^ C Ne{p) holds. If p G Homc{m, n) is H-SStable then any 
non- zero vector O < X G i?™ belongs to the (open) hyperbolic cone Ce{p)- I 

20 



We need the next fundamental fact due to L. Garding [B] (we recommend 
the very readable treatinent in 



Theorem 5.3: Let p G Homc{m,n) be e-hyperbolic polynomial and d G 
Ceip) C -R™. Then p is also d-hyperbolic and Cd{p) = Ce{p), Ndip) = Ne{p). 
Moreover cone Ce{p), called hyperbolic cone, is convex. 

Corollary 5.4: 

1. For any two vectors in the hyperbolic cone di,d2 G Ce{p) the following 
set equality holds: 

NdAp) n i-NdAp)) = NdM n {-NdM) = Nullp. (36) 

Thus Nullp C -R™ is a linear subspace. 

2. 

Nullp = {X eR^ : p{Y + X)= p{Y) for all Y G C""} (37) 

Let Pr (Nullp) be orthogonal projector on the linear subspace Nullp. It 
follows from (3^ that 

p{Y) =p{{I- Pr{Nullp))Y) (38) 

3. Let p G Hom^{n,n) be a doubly-stochastic H-Stable polynomial. If 
(?/i,---,?/n) e Nullp then 



y, + ... + y^ = 0. (39) 



Proof: 



1. It is well known and obvious that if iT is a convex cone in some hnear 
space L over reals then the intersection Kr\{—K) is a hnear subspace 
of L. 

2. Let T G Ce{p) and X G Nullp. Then aU the roots of the equation 
p{xT+X) = O are equal to zero. Since p(T) ^ O and the polynomial p G 
Homc{m,n) is homogeneous, hence p{xT + X) = x"'p{T). Therefore 
p(T + X) = p{T) for all T G Ce(jo). As Ce(p) is a non-empty open 
subset of i?"^, equality (1571) follows from the analyticity of p. 
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3. Consider the vector of all ones e = (1,...,1) G i?" and a vector Y = 
{yi, ■ ■■iVn) e Nullp. Then d{t) = p{e + tY) = p{e) for all t E R. 
Therefore 

o = d'(o)= J2 y^T^p{lA,■■■,l) = yl + ... + yn■ 

l<i<n '^^i 
I 

Example 5.5: 

1. Consider the power polynomial q G Hom^{n, n), q{xi, . . . , x„) = (aiXi+ 
. . . + ünXnY- If the non-zero vector a = (oi, . . . , a„) G i?" then the 
power polynomial q is H-Stable. The correspondind linear subspace 

Nullp = a-^ =: {(yi, . . . , |/„) G i?" : ^ OjI/j = 0}, dim{Nullp) = n-1; 

l<i<n 

and Capijp) = n^^Y[i<i<n'^i- Therefore Cap{q) 7^ O iff a G -R++- 
It is easy to see that 

Qn 

(aiXi + . . . + ünXnY = '^'öi . . .ttn- 



dxi . . . dXn 

Therefore 



-q{0,...,0) = Cap{q)-. (40) 



dxi . . . dxn ' ' n^' 

If dim{Nullp) = n — 1 and a polynomial p G Hom^ln, n) is H-Stable 
then p(a;i, . . . ,a;„) = (61X1 + . . . + bnXn)"' for some non-zero vector 
b = (bi, . . . , b„) G i?" . The power polynomial jo is doubly-stochastic iff 
6j = -, 1 < z < n. 

2. Let p G Homc{m,n) be an e-hyperbolic polynomial, D G Ce(p) C R"^ 
and X G i?™. Suppose that the univariate polynomial R{t) = p{tD + 
X) = a{t + fe)", b e R. Define the next real vector Y = -bD + X. 
Then all the roots of the equation p{Y — XD) = O are equal to zero. 
Therefore Y G Noip) R (-A^Db)) = iVw/^p- 
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5.2 Uniqueness 

Definition 5.6: We call a H-Stable polynomial p e Hom+{n,n) extremal 
if 

Cap{p) > O and 

p{0,...,0) = —Capip). (41) 



dxi . . . dxn n 

I 

Our goal is the next theorein 

Theorem 5.7: A H-Stable polynomialp G Hom^{n,n) zs extremal if and 

only if 

p{xi, . . .,Xn) = {aixi + . . . + a„x„)" 

for some positive real numbers ai, . . . , a„ > 0. 

In other words, the equality ^1\ ) holds iff dim{Nullp) = n — 1. 

Notice that the "if" part is simple and follows from the equality ^^ . 

We coUect the basic properties of extremal polynomials in the next 
proposition. 

Proposition 5.8: 

1. If a H-Stable polynomialp G IIom^{n,n) is extremal then all its 
coefficients are positive real numbers. 

2. Let c = (ci,...,c„) G i?"^ and p G IIom^{n,n). Define the scaled 
polynomial Pc as Pc{xi, . . . ,Xn) = p{ciXi, . . . , c„x„). Ifp G IIom^{n,n) 
is H-Stable extremal polynomial then also the scaled polynomial p^ 
is. 

Proof: 

1. Our goal in this step is to show that if equation (14T]) holds then the 
H-Stable polynomial p is, in fact, H-SStable and therefore has all 
positive coefficients. 

Since G{2) ■ ■ ■ G{n) = ^ and the function G is strictly decreasing on 
[0,oo), hence it follows from fITÏÏl) that degp{n) = n. Since the inequal- 
ity flT^ is invariant with respect to permutations of variables hence 
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degp{i) = n : 1 < i < n. Which means that p(ej) > O : 1 < i < n, 
where {ci, . . . , e„} is the Standard orthonorinal basis in R". Therefore 
the polynomial p is H-SStable. Thus its coefiicients are strictly posi- 
tive real numbers and all non-zero vectors Y G R^ belong to its open 
hyperbolic cone Ce{p). 

2. First, if p G Hom^{n,n) is H-Stable then clearly the scaled polyno- 
mial Pc is also H-Stable. 

It foUows from the definition of the Capacity that Cap{pc) = Ci ■ ■ ■ CnCap{p). 
We get, by a direct computation, that 

Qn Qn 

-Pc(0, . . . , 0) = Cl • • • c„- ^— P(0, . . . , 0). 



f-\ r\ lr^\ 7 3 "/ ^-L ^/t r\ r\ i 

OXi . . . OXn OXi . . . OXn 

This proves that the set of H-Stable extremal polynomials is invari- 
ant with respect to the scaling. 



We need the following simple result (it was essentially proved in Lemma 
3.8 from (8]). 

Fact 5.9: Considerp G Hom++{n,n),p{xi, . . .x„) = Eri,...,r„ ari,...,r„ I\i<i<nX?- 
Then there exists a positive vector t =: (ti, . . . , t„) G i?"^, ti ■ ■ ■ t„ = 1 such 
that 

p{ti,...,tn) =Cap{p) = ^^^inf p{xi,...,Xn). (42) 

Xi>U ,l<.l<n]Xi- ■ ■Xn = l 

Consider the corresponding scaled polynomial pt- Then the polynomial 
q = J* -^^^ is doubly-stochastic. I.e. ^^(1, 1, . . . , 1) = 1, 1 < i < n.. 

Proof: Consider a subset 

T = {(xi,. . . ,x„) G i?"+ : Xi ■ --Xn = 1 and p(xi, . . . , x„) < p(l, . . . , 1)}. 

In order to prove that the infimum is attained, it is sufficiënt to show that 
the subset T is compact. Clearly, T is closed, and we need to prove that T 
is bounded. Let (xi, . . . , x„) G T and assume WLOG that maxi<j<„ Xj = Xi. 
Then 

p(l, . . .,1) 

an,o,...,oXi < p{xi, . . . , x„) < j9(l, . . . , 1) ^ max Xj < '■ '- — < oo. 

l<ï<n a„,o,...,o 
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This shows the desired boundness of T and the existence of the mmimuiii. 
Consider a positive vector (ti, . . . , t„) such that p{ti, ■ . ■ ,tn) = Capijp). De- 
fine ai = log(tj), 1 < i < n. Then 

p{exp{ai), ..., exp{an)) = ^ min p{exp{Pi), ..., exp{Pn))- 

/5l+.../5„=0 

Therefore there exists the Lagrange multiplier 7 such that 

d d 

p{exp{ai), ..., exp{an)) = U—-p{ti, . . . ,tn) = ^,1 < i < n. 



oai öti 

It follows from the Euler's identity that 7 = p(ti, . . . , t„) and 

d d 

— g(l,...,l) = (p(ti,...,t„))"Hi— p(ti,...,t„) = 1. 



Remark 5.10: It is easy to prove that, in fact, the minimum in fH21) is 
attained uniquely. It was proved in [12] that if p G Homj^{n, n) is H-Stable 
then the minimum in (H2!) exists and attained uniquely iff 

Qn Qn 

-p(0, . . . , 0), ^ ^ ^ ^P(0, . . . , 0) > O : 1 < ï 7^ j < n. 



OXjOXj [lfnj^{i,j) ^-^171 OXiOXi [lm^(ij) OX^-, 



Proof: (Proof of Theorem dM]).) 

It follows from Proposition fIS.Sp and Fact (15. 91) that we can assume without 
loss of generality that the H-Stable extremal polynomial p G Hom^^{n, n) 
is doubly-stochastic and all its coefficients are positive real numbers. 

1. Using uniqueness part of Lemma( l3.2l) 

Let {ei, . . . , e„} be the Standard basis in i?" and p is now a H-SStable 
doubly-stochastic polynomial with positive coefficients, p satisfies the 
equality fHTl) . We need to look at the case of equality in fl35|) . Recall the 
polynomial g„_i is given by g„-i(xi, . . . , Xn-i) = ^p{xi, ..., Xn-i, 0). 

Clearly, the polynomial g„_i G Hom^^{n — l,n — l) also has positive 
coefficients. Let 

Qn-l\tl,ni ■ ■ ■ ■,tn-l,n) = mhl g„_i(Xi, . . . , X^-i) . 

xi,...,x„_i>0,ni<,<„_ia:i=l 
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The existence of such a vector was proved in Proposition (15.91) . 
It foUows from the uniqueness part of Lemma (13. 21) that the univariate 
polynomial R{t) = p(ti,n, ■ ■ ■ , tn-i,n, t) = p(X,i<i<n-i U,n-iei + ten) has 
n equal negative roots: R{t) = b{t + a„)"; a„, b > 0. This fact imphes, 
as in the second part of Example (15.51) , that 

Kn =: en - Yl ^3,n'^3 ^ Nullp, üj^n = ^^ > 0. 

Since p is doubly-stochastic, hence it foUows from (15Ü]) that the coor- 
dinates of Kn sum to zero. Which gives that J2jjtn(^j,n = 1- 

In the same way, we get that there exists an n x n column stochastic 
matrix A with the zero diagonal and the positive off-diagonal part such 
that the vectors 

Ki =: Cj — ^ ^j,i^j ^ Nullp, i < i < n. 

2. Recall that our goal is to prove that dim{Nullp) = n — 1. It follows that 

dim{Nullp) > dim{L{Ki, . . . , Kn)) = Rank{I — A), 
where L{Ki, . . . , Kn) is the minimal linear subspace containing the set 

{K,,...,Kn}. 

Since the polynomial p is non-zero thus dim{Nullp) < n — 1. It is easy 

to see that Rank{I — A) = n — 1. Indeed, any principal n — 1 x n — 1 

submatrix of J — A is strictly diagonally dominant and, therefore, is 

nonsingular. 

We finally conclude that dim,{Nullp) = n — 1. 



6 Comments 

1. Falikman [5j and Egorychev ^ publications were foliowed by a fiurry 
of expository papers, which clarified and popularized the proofs. The 
author learned the Egorychev's proof from [17j. It is our guess that 
many scientists first learned about Alexandrov inequalities for mixed 
discriminants and Alexandrov-Fenchel inequalities for mixed volumes 
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pp in one of those expository papers. We would like to distinguish the 
foUowing two papers: (Tü] and [23]. They both explicitly connected 
Alexandrov inequalities for mixed discriminants with honiogeneous hy- 
perbolic polynomials. The paper p£j was, essentially a rediscovery of 
Garding's theory [6]. Still, as the author had read [16] before reading 
[ü], the paper [Tü] gave us the first hint for the possibihty of our ap- 
proach. 

The paper [25], apparently written as a technical report in 1981 and 
pubhshed only in 2006 in an obscure book, is technically very similar 
to [Tü]. Besides, it imphcitly introduced the Bapat's conjecture. 
Other related pubhcations are D. London's (univariate) papers 



As far as we know, there were no previously pubhshed connections 
between Shrijver-Vahant conjecture, which was thought to be of purely 
combinatorial nature, and stable/hyperbohc polynomials. 

2. Two main ingredients of our approach, which make the proofs sim- 
ple, are the usage of the notion of Capacity and Lemma (13.21) . They 
together allowed the simple induction. The induction, used in this 
paper, is by partial differentiation. It is very similar to the induc- 
tive proofs of hyperbolic polynomials analogues of Alexandrov inequal- 
ities for mixed discriminants in [Iü],[23]. Using our terminology, these 
analogues correspond to the fact that the polynomial q2{xi,X2) = 
dx" dx P(^i' a;2, O, . . . , 0) is either zero or H-Stable provided the poly- 
nomial p is H-Stable. 

The idea to use Capacity in the context of permanents is implicit 
in [22j. The notion of Capacity was crucial for algorithmic results 
in [7], [8j as \og{Cap{p)) = infyj+...+y^^=olog (jo(e^i, . . . , e^")) and the 
functional log (p(e^\ . . . , e^")) is convex for any polynomial with non- 
negative coefficients. 

Probably, the papers [7] , [8] were the first to reformulate Van der Waer- 
den/Bapat conjectures as in inequality (1271) . Although quite simple, it 
happened to be a very enlighting observation. 

3. Our, inductive by the partial differentiation, approach was initiated in 
jlOj . The main tooi there was Vinnikov-Dubrovin determinantal rep- 
resentation [2Ü] of hyperbolic homogeneous polynomials in 3 variables. 
The paper [TÜ] proved the implication 
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Cap{p) > O =^ Qx %x P(Q' • • • ' 0) > O fo^ H-Stable polynomials 
p e Homj^{n,n). Additionally, it was proved that in this H-Stable 
case the functional Rankp{S) = maxa^ ,...,r„j^oJ2j£s''^j is submodular 
and 

an,...,r-„ > O ^=^ ^ r, < RankpiS) : S C {1, . . . ,n}. (43) 

The characterization ( H3l) is a far reaching generahzation of the Hall- 
Rado theorem. 

The paper p!2] provides algorithniic apphcations of these resuhs: strongly 
polynomial deterministic algorithms for the membership problem as for 
the support as well for the Newton polytope of H-Stable polynomials 
p G Hom+{m, n), given as oracles. 
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